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The status of the search for neutrinoless double beta decay is reviewed. 
The effort to reach the sensitivity needed to cover the effective Majorana 
neutrino mass corresponding to the degenerate and inverted mass hierar- 
f"^ ■ chy is described. Various issues concerning the theory (and phenomenol- 

ogy) of the relation between the 0vf3/3 decay rate and the absolute neu- 
trino mass scale are discussed, in particular the issue of mechanism of 
the Qv(3(3 decay. Finally, the relation between the neutrino magnetic 
|| moments and the charge conjugation property (Dirac vs. Majorana) is 

D . described. 



1.1. Introduction - fundamentals of (3(3 decay 

In the recent past neutrino oscillation experiments have convincingly shown 
that neutrinos have a finite mass. However, in oscillation experiments only 
the differences of squares of the neutrino masses, Am 2 = |m| — mf\, can 
be measured, and the results do not depend on the charge conjugation 
properties of neutrinos, i.e., whether they are Dirac or Majorana fermions. 
Nevertheless, a lower limit on the absolute value of the neutrino mass scale, 
m scaie — -\/|Ato 2 |, has been established in this way. Its existence, in turn, 
is causing a renaissance of enthusiasm in the double beta decay community 
which is expected to reach and even exceed, in the next generation of exper- 
iments, the sensitivity corresponding to this mass scale. Below I review the 
current status of the double beta decay and the effort devoted to reach the 
required sensitivity, as well as various issues in theory (or phenomenology) 
concerning the relation of the 0vf3(3 decay rate to the absolute neutrino 
mass scale and to the general problem of the Lepton Number Violation 
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(LNV). 

But before doing that I very briefly summarize the achievements of the 
neutrino oscillation searches and the role that the search for the neutrino- 
less double beta decay plays in the elucidation of the pattern of neutrino 
masses and mixing. In these introductory remarks I use the established 
terminology, some of which will be defined only later in the text. 

There is a consensus that the measurement of atmospheric neutrinos by 
the SuperKamiokande collaboration 1 can be only interpreted as a conse- 
quence of the nearly maximum mixing between and v r neutrinos, with 
the corresponding mass squared difference |Am„ tm | ~ 2.4 x 10~ 3 cV 2 . This 
finding was confirmed by the K2K experiment 2 that uses accelerator 
beam pointing towards the SuperKamiokande detector 250 km away, as 
well as by the very recent first result of the MINOS experiment located 
at the Sudan mine in Minnesota, 735 km away from the Fermilab. 3 Sev- 
eral large long-baseline experiments are being built to further elucidate this 
discovery, and determine the corresponding parameters more accurately. 

At the same time the "solar neutrino puzzle" , which has been with us 
for over thirty years since the pioneering chlorine experiment of Davis, also 
reached the stage where the interpretation of the measurements in terms of 
oscillations between the z/ e and some combination of the active, i.e., and 
v T neutrinos, is inescapable. In particular, the juxtaposition of the results 
of the SNO experiment 5 and SuperKamiokande, 6 together with the earlier 
solar neutrino flux determination in the gallium experiments, 7,8 leads to 
that conclusion. The value of the corresponding oscillation parameters, 
however, remained uncertain, with several "solutions" possible, although 
the so-called Large Mixing Angle (LMA) solution with sin 2 29 so i ~ 0.8 and 
Am 2 o( ~ 10~ 4 eV 2 was preferred. A decisive confirmation of the "solar" 
oscillations was provided by the nuclear reactor experiment KamLAND 9 ' 10 
that demonstrated that the flux of the reactor v e is reduced and its spectrum 
distorted at the average distance ~ 180 km from nuclear reactors. 

The pattern of neutrino mixing is further simplified by the constraint 
due to the Chooz and Palo Verde reactor neutrino experiments 11,12 which 
lead to the conclusion that the third mixing angle, 6*13, is small, sin 2 2#i3 < 
0.1. The two remaining possible neutrino mass patterns are illustrated in 
FigO 

Altogether, clearly a lower limit for at least one of the neutrino masses, 
y/ Am„ tm ~ 0.05 eV has been established. However, the oscillation ex- 
periments cannot determine the absolute magnitude of the masses and, in 
particular, cannot at this stage separate two rather different scenarios, the 
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Fig. 1.1. Schematic illustration (mass intervals not to scale) of the decomposition of 
the neutrino mass eigenstates in terms of the flavor eigenstates. The two hierarchies 
cannot be, at this time, distinguished. The small admixture of v e into v$ is an upper 
limit. 



hierarchical pattern of neutrino masses in which m ~ v Am 2 and the de- 
generate pattern in which m 3> V Am 2 . It is hoped that the search for 
the neutrinoless double beta decay, reviewed here, will help in foreseeable 
future in determining, or at least narrowing down, the absolute neutrino 
mass scale, and in deciding which of these two possibilities is applicable. 

Moreover, the oscillation results do not tell us anything about the prop- 
erties of neutrinos under charge conjugation. While the charged leptons 
are Dirac particles, distinct from their antiparticles, neutrinos may be the 
ultimate neutral particles, as envisioned by Majorana, that are identical to 
their antiparticles. That fundamental distinction becomes important only 
for massive particles and becomes irrelevant in the massless limit. Neutrino- 
less double beta decay proceeds only when neutrinos are massive Majorana 
particles, hence its observation would resolve the question. 

Double beta decay (f3f3) is a nuclear transition (Z,A) — > (Z + 2, A) in 
which two neutrons bound in a nucleus are simultaneously transformed into 
two protons plus two electrons (and possibly other light neutral particles). 
This transition is possible and potentially observable because nuclei with 
even Z and N are more bound than the odd-odd nuclei with the same 
A = N + Z. Analogous transition of two protons into two neutrons are 
also, in principle, possible in several nuclei, but phase space considerations 
give preference to the former mode. 

An example is shown in Fig. 11.21 The situation shown there is not 
exceptional. There are eleven analogous cases (candidate nuclei) with the 
Q- value (i.e., the energy available to leptons) in excess of 2 MeV. 
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There are two basic modes of the (3(3 decay. In the two-neutrino mode 
(2v(3(3) there are two v e emitted together with the two e~. Lepton number 
is conserved and this mode is allowed in the standard model of electroweak 
interaction. It has been repeatedly observed by now in a number of cases 
and proceeds with a typical half- life of ~ 10 20 years. In contrast, in the 
neutrinoless mode (0^/3/3) only the 2e~ are emitted and nothing else. That 
mode clearly violates the law of lepton number conservation and is forbid- 
ded in the standard model. Hence, its observation would be a signal of a 
"new physics". 
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Fig. 1.2. Atomic masses of the isotopes with A = 136. Nuclei 136 Xe, 136 Ba and 136 Cc 
are stable against the ordinary fi decay; hence they exist in nature. However, energy 
conservation alone allows the transition 136 Xc — > 136 Ba + 2e~ (+ possibly other neutral 
light particles) and the analogous decay of 136 Ce with the positron emission. 

The two modes of the /3/3 decay have some common and some distinct 
features. The common features are: 

• The leptons carry essentially all available energy. The nuclear recoil 
is negligible, Q/Am p <C 1. 

• The transition involves the + ground state of the initial nucleus 
and (in almost all cases) the + ground state of the final nucleus. In 
few cases the transition to an excited + state in the final nucleus 
is energetically possible, but suppressed by the smaller phase space 
available. (But the 2vf3j3 decay to the excited + state has been 
observed in few cases.) 
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Both processes are of second order of weak interactions, ~ G F , 
hence inherently slow. The phase space consideration alone (for the 
2v/3(3 mode ~ Q 11 and for the 0v/3(3 mode ~ Q 5 ) give preference 
to the Qvj3f3 which is, however, forbidden by the lepton number 
conservation. 



The distinct features are: 

• In the 2v(3j3 mode the two neutrons undergoing the transition are 
uncorrelated (but decay simultaneously) while in the §v(3(3 the two 
neutrons are correlated. 

• In the 2vj3[3 mode the sum electron kinetic energy Ti+T% spectrum 
is continuous and peaked below Q/2. As T1 + T2 — ► Q the spectrum 
approaches zero approximately like (AE/Q) 6 . 

• On the other hand, in the 0^/3/3 mode T\ + T2 = Q smeared only 
by the detector resolution. 

These last features allow one to separate the two modes experimentally 
by observing the sum electron spectrum with a good energy resolution, even 
if the corresponding decay rate for the 0v(3j3 mode is much smaller than for 
the 2vf3j3 mode. This is illustrated in Fig |1.3l where the insert that includes 
the Of peak and the 2v tail shows the situation for the rate ratio of 1 : 10 6 
corresponding to the most sensitive current experiments. 
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Fig. 1.3. Separating the 0vf3f3 mode from the 2f/3/3 by the shape of the sum elec- 
tron spectrum (kinetic energy K e of the two electrons), including the effect of the 2% 
resolution smearing. The assumed 2u/0u rate ratio is 10 2 , and 10 6 in the insert. 
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Various aspects, both theoretical and experimental, of the (3(3 decay 
have been reviewed many times. Here I quote just the more recent review 
articles, 13-16 earlier references can be found there. 

In this introductory section let me make only few general remarks. The 
existence of the Qv(3(3 decay would mean that on the elementary particle 
level a six fermion lepton number violating amplitude transforming two u 
quarks into two d quarks and two electrons is nonvanishing. As was first 
pointed out by Schechter and Valle 17 more than twenty years ago, this fact 
alone would guarantee that neutrinos arc massive Majorana fermions (see 
Fig ll.4p . This qualitative statement (or theorem), unfortunately, does not 
allow us to deduce the magnitude of the neutrino mass once the rate of the 
Qv(3(3 decay have been determined. 




Fig. 1.4. By adding loops involving only standard weak interaction processes the 0vf3f3 
decay amplitude (the black box) implies the existence of the Majorana neutrino mass. 

There is no indication at the present time that neutrinos have nonstan- 
dard interactions, i.e. they seem to have only interactions carried by the 
W and Z bosons that are contained in the Standard Electroweak Model. 
All observed oscillation phenomena can be understood if one assumes that 
that neutrinos interact exactly the way the Standard Model prescribes, but 
are massive fermions forcing a generalization of the model. If we accept 
this, but in addition assume that neutrinos are Majorana particles, we can 
in fact relate the Qv(3(3 decay rate to the quantity related to the absolute 
neutrino mass. With these caveats that relation can be expressed as 

-±- = G o »(Q,Z)\M°r{m fj0 r , (1.1) 

1 l/2 

where G 0u (Q,Z) is a phase space factor that depends on the transition 
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Q value and through the Coulomb effect on the emitted electrons on the 
nuclear charge Z and that can be easily and accurately calculated, M 0l/ 
is the nuclear matrix element that can be evaluated in principle, although 
with a considerable uncertainty, and finally the quantity (mpp) is the effec- 
tive neutrino Majorana mass, representing the important particle physics 
ingredient of the process. 

In turn, the effective mass (mpp) is related to the mixing angles Oij 
that are determined or constrained by the oscillation experiments, to the 
absolute neutrino masses m, of the mass eigenstates Vi and to the as of 
now totally unknown additional parameters as fundamental as the mixing 
angles 9ij, the so-called Majorana phase a(i), 

(m w ) = \Ei\U ei \ 2 e ia ^mi\ . (1.2) 

Here U e i are the matrix elements of the first row of the neutrino mixing 
matrix. 

It is straightforward to use the eq. (|1.2|) and the known neutrino oscil- 
lation results in order to relate (mpp) to other neutrino mass dependent 
quantities. This is illustrated in Fig ll.51 Traditionally such plot is made 
as in the left panel. However, the lightest neutrino mass m m i n is not an 
observable quantity. For that reason the other two panels show the relation 
of (mpp) to the sum of the neutrino masses M = Sm, and also to (mp) 
that represents the parameter that can be determined or constrained in 
ordinary (3 decay, 

(m fj ) 2 =S i |f/ ei | 2 m? • (1-3) 

Several remarks are in order. First, the observation of the 01/0(3 decay 
and determination of (mpp), even when combined with the knowledge of M 
and/or (mp) does not allow, in general, to distinguish between the normal 
and inverted mass orderings. This is a consequence of the fact that the 
Majorana phases are unknown. In regions in Fig. 1 1.51 where the two hatched 
bands overlap it is clear that two solutions with the same (mpp) and the 
same M (or the same (mp)) exist and cannot be distinguished. 

On the other hand, obviously, if one can determine that (mpp) > 0.1 
eV we would conclude that the mass pattern is degenerate. And in the so 
far hypothetical case that one could show that (mpp) < 0.01 - 0.02 eV but 
nonvanishing nevertheless the normal hierarchy would be established. 

a In that case also the (mp) in the right panel would not represent the quatity directly 
related to the ordinary (3 decay. There are no realistic ideas, however, how to reach the 
corresponding sensitivity in ordinary (3 decay at this time. 
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It is worthwhile noting that if the inverted mass ordering is realized 
in nature, (and neutrinos are Majorana particles) the quantity (mpp} is 
constrained from below by ~ 0.01 eV. This value is within reach of the 
next generation of experiments. Also, in principle, in the case of the normal 
hierarchy while all neutrinos could be massive Majorana particles it is still 
possible that {mpp) = 0. Such a situation, however, requires "fine tuning" 
or reflects a symmetry of some kind. 




-3 -? -1 -7 -1 -3 -? -1 

10 10 10 10 10 10 10 10 

m m!n / eV M / eV <m„>/eV 



Fig. 1.5. The left panel shows the dependence of (mpp) on the mass of the lightest neu- 
trino m m i„, the middle one shows the relation between (mpp) and the sum of neutrino 
masses M = T,mi determined or constrained by the "observational cosmology" , and the 
right one depicts the relation between (mpp) and the effective mass (mp) determined 
or constrained by the ordinary j3 decay. In all panels the width of the hatched area is 
due to the unknown Majorana phases and therefore irreducible. The solid lines indicate 
the allowed regions by taking into account the current uncertainties in the oscillation 
parameters; they will shrink as the accuracy improves. The two sets of curves correspond 
to the normal and inverted hierarchies, they merge above about (mpp) > 0.1 eV, where 
the degenerate mass pattern begins. 



Let us remark that the 0vf3(3 decay is not the only LNV process for 
which important experimental constraints exist. Examples of the other 
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analogous processes are 

yT + (Z, A) — > e+ + (Z - 2, A); exp. branching ratio < lCT 12 , 

K + — ► /i + /i + 7r~; exp. branching ratio < 3 x 10 -9 , 

u e emission from the Sun; exp. branching ratio < 10~ 4 . (1-4) 

However, detailed analysis suggests that the study of the 0v(3(3 decay is 
by far the most sensitive test of LNV. In simple terms, this is caused by 
the amount of tries one can make. A 100 kg Of/3/3 decay source contains 
~ 10 27 nuclei. This can be contrasted with the possibilities of first pro- 
ducing muons or kaons, and then searching for the unusual decay channels. 
The Fermilab accelerators, for example, produce "a few" x 10 20 protons on 
target per year in their beams and thus correspondingly smaller numbers 
of muons or kaons. 

1.2. Mechanism of the 0v/3f3 decay 

It has been recognized long time ago that the relation between the Oi//30 de- 
cay rate and the effective Majorana mass {nipp) is to some extent problem- 
atic. The assumption leading to the eq. (|l.l|) is rather conservative, namely 
that there is an exchange of a virtual light, but massive, Majorana neu- 
trino between the two nucleons undergoing the transition, and that these 
neutrinos interact by the standard left-handed weak currents. However, 
that is not the only possible mechanism. LNV interactions involving so far 
unobserved heavy (~ TeV) particles can lead to a comparable 0^/3/3 decay 
rate. Thus, in the absence of additional information about the mechanism 
responsible for the 0^/3/3 decay, one could not unambiguously infer (rapp) 
from the 0^/3/3 decay rate. 

In general 0^/3/3 decay can be generated by (i) light massive Majorana 
neutrino exchange or (ii) heavy particle exchange (see, e.g. Refs. 18 ' 19 ), re- 
sulting from LNV dynamics at some scale A above the electroweak one. 
The relative size of heavy (Ah) versus light particle (A£) exchange contri- 
butions to the decay amplitude can be crudely estimated as follows: 20 

A r i (mgg) , r 2 A H M^(fc 2 ) 

where {mpp) is the effective neutrino Majorana mass, (fc 2 ) ~ (50 MeV) 2 
is the typical light neutrino virtuality, and A is the heavy scale relevant to 
the LNV dynamics. Therefore, A H /A L - 0(1) for {m pfi ) - 0.1 - 0.5 eV 
and A ~ 1 TeV, and thus the LNV dynamics at the TeV scale leads to 
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similar Of/3/3 decay rate as the exchange of light Majorana neutrinos with 
the effective mass (TO/3/3) ~ O.f — 0.5 eV. 

Obviously, the lifetime measurement by itself does not provide the 
means for determining the underlying mechanism. The spin-flip and non- 
flip exchange can be, in principle, distinguished by the measurement of the 
single-electron spectra or polarization (see e.g., 21 ). However, in most cases 
the mechanism of light Majorana neutrino exchange, and of heavy particle 
exchange cannot be separated by the observation of the emitted electrons. 
Thus one must look for other phenomenological consequences of the differ- 
ent mechanisms other than observables directly associated with 0v(3j3. Here 
1 discuss the suggestion 22 that under natural assumptions the presence of 
low scale LNV interactions also affects muon lepton flavor violating (LFV) 
processes, and in particular enhances the fi — > e conversion compared to 
the ji — ► e7 decay. 

The discussion is concerned mainly with the branching ratios B^ e -y = 
r(/i — > e^)/T^ and B^ e = r conv /r capt , where [i — > is normalized 
to the standard muon decay rate r^ "* = (G 2 ? to 5 J/(1927t 3 ), while \i — ► e 
conversion is normalized to the corresponding capture rate r cap t. The main 
diagnostic tool in the analysis is the ratio 



and the relevance of our observation relies on the potential for LFV discov- 
ery in the forthcoming experiments MEG 23 Qi -> ey) and MECO 24 (fi -> e 
conversion^ that plan to improve the current limits by several orders of 
magnitude. 

It is useful to formulate the problem in terms of effective low energy 
interactions obtained after integrating out the heavy degrees of freedom 
that induce LNV and LFV dynamics. If the scales for both LNV and LFV 
are well above the weak scale, then one would not expect to observe any 
signal in the forthcoming LFV experiments, nor would the effects of heavy 
particle exchange enter 0v(3f3 at an appreciable level. In this case, the only 
origin of a signal in 0v(3f3 at the level of prospective experimental sensitivity 
would be the exchange of a light Majorana neutrino, leading to eq. (|l.l|L 
and allowing one to extract (TO/3/3) from the decay rate. 

In general, however, the two scales may be distinct, as in SUSY-GUT 25 
or SUSY see-saw 26 models. In these scenarios, both the Majorana neutrino 
mass as well as LFV effects are generated at the GUT scale. The effects of 

b Even though MECO experiment was recently cancelled, proposals for experiments with 
similar sensitivity exist elsewhere. 



— B^ e /B l 



fjt,— ye'y ) 



(1.6) 
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heavy Majorana neutrino exchange in Qv(3f3 are, thus, highly suppressed. In 
contrast, the effects of GUT-scale LFV are transmitted to the TeV-scale by 
a soft SUSY-breaking sector without mass suppression via renormalization 
group running of the high-scale LFV couplings. Consequently, such scenar- 
ios could lead to observable effects in the upcoming LFV experiments but 
with an 0(a) suppression of the branching ratio B^ e relative to B^ ei 
due to the exchange of a virtual photon in the conversion process rather 
than the emission of a real one, thus 1Z <~ 10~( 2-3 ) in this case. 

The case where the scales of LNV and LFV are both relatively low (~ 
TeV) is more subtle. This is the scenario which might lead to observable 
signals in LFV searches and at the same time generate ambiguities in in- 
terpreting a positive signal in 0v/3[3. Therefore, this is the case where one 
needs to develop some discriminating criteria. 

Denoting the new physics scale by A, one has a LNV effective lagrangian 
of the form 

Coufffi = Yl % & & = FiQ ^ 2 qeT 3 e c , (1.7) 

i 

where we have suppressed the flavor and Dirac structures (a complete list 
of the dimension nine operators Oi can be found in Ref. 19 ). 
For the LFV interactions, one has 

£ LFV = £ || Oi, (1.8) 

i 

and a complete operator basis can be found in Rcfs. 27,28 The LFV operators 
relevant to our analysis are of the following type (along with their analogues 
with L <-» R): 

O aL = j^I^a^ipljL F» v + h.c. 
Oil = T~Ll c jL Wl^l 

O eq = JlTtij qT q q . (1.9) 

Operators of the type O a are typically generated at one- loop level, hence 
our choice to explicitly display the loop factor 1/(4tt) 2 . On the other hand, 
in a large class of models, operators of the type Oe or Oe q are generated by 
tree level exchange of heavy degrees of freedom. With the above choices, 
all non-zero Ci and Ci are nominally of the same size, typically the prod- 
uct of two Yukawa-like couplings or gauge couplings (times flavor mixing 
matrices). 
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With the notation established above, the ratio TZ of the branching ratios 
/i — » e to [i — > e+7 can be written schematically as follows (neglecting flavor 
indices in the effective couplings and the term with L <-> R): 



TZ 



A 2 

Ai e 2 c aL + e 2 (A 2 C£ L + A 3 c f? ) log 



TO 2 



48tt 2 

+ A 4 (47r) 2 c, g + ..."/ [e 2 (|c ffL | 2 + | C(Tfl | 2 )] . (1.10) 

In the above formula Ai 2,3,4 are numerical factors of 0(1), while the overall 
factor arises from phase space and overlap integrals of electron and 
muon wavefunctions in the nuclear field. For light nuclei <& = ( Z F 2 ) / (g v + 
3g\) ~ 0(1) (gv,A are the vector and axial nucleon form factors at zero 
momentum transfer, while F p is the nuclear form factor at q 2 — — m 228 ). 
The dots indicate subleading terms, not relevant for our discussion, such as 
loop- induced contributions to ct and ci q that are analytic in external masses 
and momenta. In contrast the logarithmically-enhanced loop contribution 
given by the second term in the numerator of 1Z plays an essential role. This 
term arises whenever the operators Oil,r and/or Og q appear at tree-level 
in the effective theory and generate one-loop renormalization of Oi q 27 (see 
Fig. Q3J . 



O e ,0 




f = t,q 



Fig. 1.6. Loop contributions to fi — * e conversion through insertion of operators Oi or 
Oi q , generating the large logarithm. 

The ingredients in eq. (jl.lOp lead to several observations: (i) In absence 
of tree-level cm, and a q , one obtains 1Z ~ ($ A 2 a)/(12w) ~ 10 -3 — 10 -2 , 
due to gauge coupling and phase space suppression, (ii) When present, the 
logarithmically enhanced contributions compensate for the gauge coupling 
and phase space suppression, leading to 1Z ~ 0(1). (ifi) If present, the 
tree-level coupling cg q dominates the fj, — > e rate leading to TZ S> 1. 

Thus, we can formulate our main conclusions regarding the discriminat- 
ing power of the ratio TZ: 
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(1) Observation of both the LFV muon processes fi — > e and fj, — ► &y 
with relative ratio 1Z ~ 10~ 2 implies, under generic conditions, that 
To v pp ~ {mpp) 2 ■ Hence the relation of the 0v(3j3 lifetime to the absolute 
neutrino mass scale is straightforward. 

(2) On the other hand, observation of LFV muon processes with relative 
ratio 1Z 3> 1CP 2 could signal non-trivial LNV dynamics at the TeV 
scale, whose effect on Ois/3/3 has to be analyzed on a case by case basis. 
Therefore, in this scenario no definite conclusion can be drawn based 
on LFV rates. 

(3) Non-observation of LFV in muon processes in forthcoming experiments 
would imply either that the scale of non-trivial LFV and LNV is above 
a few TeV, and thus T^pp ~ (mpp) 2 , or that any TeV-scale LNV is 
approximately flavor diagonal. 

The above statements are illustrated using two explicit cases: 22 the 
minimal supersymmetric standard model (MSSM) with R-parity violation 
(RPV-SUSY) and the Left-Right Symmetric Model (LRSM). 

RPV SUSY — If one does not impose R-parity conservation [R = 
(-l) 3(s ~ L)+2s ], the MSSM superpotential includes, in addition to the stan- 
dard Yukawa terms, lepton and baryon number violating interactions, com- 
pactly written as (see e.g., 29 ) 

W RPV = KjkLiL-jEl + X i]k LiQjD c k + X'l^U^D^Dl 

+^L l H u , (1.11) 

where L and Q represent lepton and quark doublet superfields, while E c , 
U c , D c are lepton and quark singlet superfields. The simultaneous pres- 
ence of A' and A" couplings would lead to an unacceptably large proton 
decay rate (for SUSY mass scale Asusy ~ TeV), so we focus on the case 
of A" = and set fi' = without loss of generality. In such case, lepton 
number is violated by the remaining terms in Wrpv, leading to short dis- 
tance contributions to Qvj3(3 [e.g., Fig J1.7f a)]. with typical coefficients [cf. 
eq. (TOD] 

Si na s Ai 2 n _ Tra 2 Ai 2 n 

A R ~ 4 ' 4 ' l 1 - 1Z J 

A J rrig m~ m x m~ 

where a s ,ot2 represent the strong and weak gauge coupling constants, re- 
spectively. The RPV interactions also lead to lepton number conserving 
but lepton flavor violating operators [e.g. Fig. I1.7f b)]. with coefficients [cf. 
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eq. (HU)] 
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AA* A'A'* 

,2 ' 



iri- 



(1.13) 



where the flavor combinations contributing to Co- can be found in Ref. 30 
Hence, for generic flavor structure of the couplings A and A' the underlying 
LNV dynamics generate both short distance contributions to Qv(3(5 and 
LFV contributions that lead to TZ > 10~ 2 . 

Existing limits on rare processes strongly constrain combinations of 
RPV couplings, assuming Asusy is between a few hundred GeV and ~ 
1 TeV. Non-observation of LFV at future experiments MEG and MECO 
could be attributed either to a larger Asusy (> few TeV) or to suppression 
of couplings that involve mixing among first and second generations. In the 
former scenario, the short distance contribution to §v(3(3 does not compete 
with the long distance one [see eq. (|1.5|) ]. so that Tqu/3/3 ~ {mpp) 2 - On the 
other hand, there is an exception to this " diagnostic tool" . If the A and A' 
matrices are nearly flavor diagonal, the exchange of superpartners may still 
make non-negligible contributions to Qu(3[3 without enhancing the ratio TZ 
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Fig. 1.7. Gluino exchange contribution to 0^/3/3 (a), and typical tree-level contribution 
to O iq (b) in RPV SUSY. 



LRSM — The LRSM provides a natural scenario for introducing non- 
sterile, right-handed neutrinos and Majorana masses. 31 The corresponding 
electroweak gauge group SU{2)l x SU(2)r x U(1)b-l, breaks down to 
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SU{2) L x U(1) Y at the scale A > O(TeV). The symmetry breaking is 
implemented through an extended Higgs sector, containing a bi-doublet $ 
and two triplets R , whose leptonic couplings generate both Majorana 
neutrino masses and LFV involving charged leptons: 



£^ pt = - L 



' (yg* + (1-14) 
- (Hp y%K L L{ - (Z^F y%A R L R . 

Here <!> = <J2$*G2, &L,r — i<^2^L,R, and leptons belong to two isospin 
doublets L l L R = [v % L Rl t % L R ). The gauge symmetry is broken through the 
VEVs (A° R ) = vr, (A° L ) = 0, ($) = diag(Ki,K 2 ). After diagonalization of 
the lepton mass matrices, LFV arises from both non-diagonal gauge inter- 
actions and the Higgs Yukawa couplings. In particular, the A^^-lepton 
interactions are not suppressed by lepton masses and have the structure 
C ~ A^\l^hij (1 ± J5)ij + h.c. The couplings are in general non- 
diagonal and related to the heavy neutrino mixing matrix. 32 

Short distance contributions to Qvj3f3 arise from the exchange of both 
heavy us and (Fig ll.Sa ). with 



52 1 92 



A5 M^ R M VR ' M^ r Ml 



(1.15) 



where is the weak gauge coupling. LFV operators are also generated 
through non-diagonal gauge and Higgs vertices, with 32 (Fig ll.8b ) 



A 2 ~ ~m^~ A 2 ~ ~MY 



e,fj,,r . (1-16) 



Note that the Yukawa interactions needed for the Majorana neutrino mass 
necessarily imply the presence of LNV and LFV couplings hij and the cor- 
responding LFV operator coefficients eg, leading to 1Z ~ 0(1)- Again, 
non-observation of LFV in the next generation of experiments would typ- 
ically push A into the multi-TeV range, thus implying a negligible short 
distance contribution to 0v/3f3. As with RPV-SUSY, this conclusion can be 
evaded by assuming a specific flavor structure, namely yu approximately 
diagonal or a nearly degenerate heavy neutrino spectrum. 

In both of these phenomenologically viable models that incorporate 
LNV and LFV at low scale (~ TeV), one finds TZ > 1CT 2 . 27 . 30 < 32 It is likely 
that the basic mechanism at work in these illustrative cases is generic: low 
scale LNV interactions (AL = ±1 and/or AL = ±2), which in general con- 
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Fig. 1.8. Typical doubly charged Higgs contribution to 0u(3/3 (a) and to Oi (b) in the 
LRSM. 

tribute to Oz//3/3, also generate sizable contributions to fi — > e conversion, 
thus enhancing this process over fi — > ej. 

In conclusion, the above considerations suggest that the ratio 1Z = 
B^e/ B^ ei of muon LFV processes will provide important insight about 
the mechanism of neutrinoless double beta decay and the use of this process 
to determine the absolute scale of neutrino mass. Assuming observation of 
LFV processes in forthcoming experiments, if 72. ~ 1CP 2 the mechanism of 
Qv(3(3 is light Majorana neutrino exchange; if 7Z 3> 10~ 2 , there might be 
TeV scale LNV dynamics, and no definite conclusion on the mechanism of 
Qv(3(3 can be drawn based only on LFV processes. 

1.3. Overview of the experimental status of search for /3{3 



The field has a venerable history. The rate of the 2v(3f3 decay was first 
estimated by Maria Goeppert-Mayer already in 1937 in her thesis work 
suggested by E. Wigner, basically correctly. Yet, first experimental obser- 
vation in a laboratory experiment was achieved only in 1987, fifty years 
later. Why it took so long? As pointed out above, the typical half-life of 
the 2vf3j3 decay is ~ 10 20 years. Yet, its "signature" is very similar to nat- 
ural radioactivity, present to some extent everywhere, and governed by the 
half-life of ~ 10 10 years. So, background suppression is the main problem 
to overcome when one wants to study either of the (3(3 decay modes. 

During the last two decades the 2v(3(3 decay has been observed in "live" 
laboratory experiments in many nuclei, often by different groups and using 
different methods. That shows not only the ingenuity of the experimen- 
talists who were able to overcome the background nemesis, but makes it 
possible at the same time to extract the corresponding 2v nuclear matrix 
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element from the measured decay rate. In the 2v mode the half-life is given 

by 

1/T 1/2 = G 2 »(Q,Z)\M 2 "\ 2 , (1.17) 

where G 2v (Q, Z) is an easily and accurately calculable phase space factor. 

The resulting nuclear matrix elements M 2v ', which have the dimension 
energy -1 , are plotted in Fig ll.9l Note the pronounced shell dependence; 
the matrix element for 100 Mo is almost ten times larger than the one for 
130 Te. Evaluation of these matrix elements, to be discussed below, is an 
important test for the nuclear theory models that aim at the determination 
of the analogous but different quantities for the §v neutrinoless mode. 




Fig. 1.9. Nuclear matrix elements for the 2v/3/3 decay extracted from the measured 
half-lives. 

The challenge of detecting the 0v(3f3 decay is, at first blush, easier. 
Unlike the continuous 2u/3/3 decay spectrum with a broad maximum at 
rather low energy where the background suppression is harder, the Ov/3/3 
decay spectrum is sharply peaked at the known Q value (see Fig |1.3[) . at 
energies that are not immune to the background, but a bit less difficult to 
manage. However, as also indicated in Fig |1.31 to obtain interesting results 
at the present time means to reach sensitivity to the Ov half-lives that are 
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~ 10 6 times longer than the 2v decay half-life of the same nucleus. So 
the requirements of background suppression are correspondingly even more 
severe. 

The historical lessons are illustrated in Fig ll.lOl where the past limits on 
the 0vf3(3 decay half-lives of various candidate nuclei are translated using 
the eq. ljl.ip into the limits on the effective mass (mpp). When plotted in 
the semi-log plot this figure represents the "Moore's law" of double beta 
decay, and indicates that, provided that the past trend will continue, the 
mass scale corresponding to Am 2 tm will be reached in about 10 years. 
This is also the time scale of significant experiments these days. Indeed, 
as discussed further, preparations are on the way to reach this sensitivity 
goal. Note that the figure was made using some assumed values of the 
corresponding nuclear matrix elements, without including their uncertainty. 
For such illustrative purposes they are, naturally, irrelevant. 

The past search for the neutrinoless double beta decay, illustrated in 
Fig ll.101 was driven by the then current technology and the resources of 
the individual experiments. The goal has been simply to reach sensitivity 
to longer and longer half-lives. The situation is different, however, now. 
The experimentalists at the present time can and do use the knowledge 
summarized in Fig Jl .51 to gauge the aim of their proposals. Based on that 
figure, the range of the mass parameter (nipp) can be divided into three 
regions of interest. 

• The degenerate mass region where all m, 3> y Am 2 tm . In that region 
{inpp) > 0.1 eV, corresponding crudely to the Qv half-lives of io 26 ~ 27 
years. To explore it (in a realistic time frame), ~ 100 kg of the decay- 
ing nucleus is needed. Several experiments aiming at such sensitivity 
are being built and should run and give results within the next 3-5 
years. Moreover, this mass region (or a substantial part of it) will be 
explored, in a similar time frame, by the study of ordinary (3 decay 
(in particular of tritium) and by the observational cosmology. These 
techniques are independent on the Majorana nature of neutrinos. It is 
easy, but perhaps premature, to envision various scenarios depending 
on the possible outcome of these measurements. 

• The so-called inverted hierarchy mass region where 20 < {mpp) < 100 
mcV and the Qv(3(3 half-lives are about 10 27 ~ 28 years. (The name is 
to some extent a misnomer. In that interval one could encounter not 
only the inverted hierarchy but also a quasi-degenerate but normal 
neutrino mass ordering. Successful observation of the 0i//3/3 decay will 
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1940 1960 1980 2000 2020 
Year 

Fig. 1.10. The limit of the effective mass (mpp) extracted from the experimental lower 
limits on the Ov/3/3 decay half-life versus the corresponding year. The gray band near 

bottom indicates the y Am^ tm value. Figure originally made by S. Elliott. 

not be able to distinguish these possibilities, as I argued above. This 
is so not only due to the anticipated experimental accuracy, but more 
fundamentally due to the unknown Majorana phases.) To explore this 
mass region, ~ ton size sources would be required. Proposals for the 
corresponding experiments exist, but none has been funded as yet, and 
presumably the real work will begin depending on the experience with 
the various ~ 100 kg size sources. Timeline for exploring this mass 
region is ~ 10 years. 
• Normal mass hierarchy region where (mpp) < 10-20 meV. To explore 
this mass region, ~ 100 ton sources would be required. There are no 
realistic proposals for experiments of this size at present. 

Over the last two decades, the methodology for double beta decay ex- 
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pcriments has improved considerably. Larger amounts of high-purity en- 
riched parent isotopes, combined with careful selection of all surrounding 
materials and using deep-underground sites have lowered backgrounds and 
increased sensitivity. The most sensitive experiments to date use 76 Ge, 
100 Mo, 116 Cd, 130 Te, and 136 Xe. For 76 Ge the lifetime limit reached im- 
pressive values exceeding 10 25 years. 33,34 The experimental lifetime limits 
have been interpreted to yield effective neutrino mass limits typically a few 
eV and in 76 Ge as low as 0.3 - 1.0 eV (the spread reflects an estimate of the 
uncertainty in the nuclear matrix elements). The sum electron spectrum 
obtained in the Heidelberg-Moscow 33 experiment is shown in Fig ll.lll over 
a broad energy range, and in Fig ll. 121 over a narrower range in the vicinity 
of the 0v Q value of 2039 keV. Some residual natural radioactivity back- 
ground lines are clearly visible in both figures, and no obvious peak at the 
0v expected position can be seen in Fig ] 1.1 21 



■l 




energy [keV] 



Fig. 1.11. The spectrum recorded in the Heidelberg-Moscow /3/3 decay experiment on 
76 Ge. Identified 7 lines are indicated. 
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Fig. 1.12. Spectrum of the Heidelberg-Moscow experiment in the vicinity of the Ou/3/3 
decay value of 2039 keV. 



Nevertheless, a subset of members of the Heidelberg-Moscow collab- 
oration reanalyzed the data (and used additional information, e.g. the 
pulse-shape analysis and a different algorithm in the peak search) and 
claimed to observe a positive signal corresponding to the effective mass of 
(mpp) = 0.391q 28 e V. 35 That report has been followed by a lively discus- 
sion. Clearly, such an extraordinary claim with its profound implications, 
requires extraordinary evidence. It is fair to say that a confirmation, both 
for the same 76 Ge parent nucleus, and better yet also in another nucleus 
with a different Q value, would be required for a consensus. In any case, 
if that claim is eventually confirmed, the degenerate mass scenario will be 
implicated, and eventual positive signal in the analysis of the tritium f3 
decay and/or observational cosmology should be forthcoming. For the neu- 
trinoless (3[3 decay the next generation of experiments, which use ~ 100 kg 
of decaying isotopes will, among other things, test this recent claim. 

It is beyond the scope of these lecture notes to describe in detail the 
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forthcoming Qv(3(3 decay experiments. Rather detailed discussion of them 
can be found e.g. in Ref. 16 Also, the corresponding chapter of the APS 
neutrino study 36 has various details. Nevertheless, let me briefly comment 
on the most advanced of the forthcoming ~ 100 kg source experiments 
CUORE, GERDA, EXO, and MAJORANA. All of them are designed to ex- 
plore all (or at least most) of the degenerate neutrino mass region {mpp) > 
0.1 eV. If their projected efficiencies and background projections are con- 
firmed, all of them plan to consider scaling up the decaying mass to ~ ton 
and extend their sensitivity to the "inverted hierarchy" region. 

These experiments use different nuclei as a source, 76 Ge for GERDA and 
MAJORANA, 130 Tc for CUORE, and 136 Xe for EXO. The requirement of 
radiopurity of the source material and surrounding auxiliary equipment is 
common to all of them, as is the placement of the experiment deep under- 
ground to shield against cosmic rays. The way the electrons are detected 
is, however, different. While the germanium detectors with their superb 
energy resolution have been used for the search of the Oi//3/3 decay for a 
long time, the cryogenic detectors in CUORE use the temperature increase 
associated with an event in the very cold TeC>2 crystals, and in the EXO 
experiment a Time Projection Chamber (TPC) uses both scintillation and 
ionization to detect the events. The EXO experiment in its final form (still 
under development and very challenging) would use a positive identification 
of the final Ba + ion as an ultimate background rejection tool. These four 
experiments are in various stages of funding and staging. First results are 
expected in about 3 years, and substantial results within 3-5 years in all of 
them. 

1.4. Nuclear matrix elements 

It follows from eq. (|1 . 1|) that (i) values of the nuclear matrix elements M 0v 
are needed in order to extract the effective neutrino mass from the measured 
Ov/3/3 decay rate, and (ii) any uncertainty in Al 0l/ causes a corresponding 
and equally large uncertainty in the extracted {mpp) value. Thus, the issue 
of an accurate evaluation of the nuclear matrix elements attracts consider- 
able attention. 

To see qualitatively where the problems are, let us consider the so- 
called closure approximation, i.e. a description in which the second order 
perturbation expression is approximated as 



M Qv = <* /mai |O (0i/) |f. 



initial 



(1.18) 
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Now, the challenge is to use an appropriate many-body nuclear model 
to describe accurately the wave functions of the ground states of the initial 
and final nuclei, {^initial) and \^ final), as well as the appropriate form of 
the effective transition operator O^ 0v ^ that describes the transformation of 
two neutrons into two protons correlated by the neutrino propagator, and 
consistent with the approximations inherent to the nuclear model used. 

Common to all methods is the description of the nucleus as a system of 
nucleons bound in the mean field and interacting by an effective residual 
interaction. The used methods differ as to the number of nucleon orbits (or 
shells and subshells) included in the calculations and the complexity of the 
configurations of the nucleons in these orbits. The two basic approaches 
used so far for the evaluation of the nuclear matrix elements for both the 2v 
and Ov 0/3 decay modes are the Quasiparticle Random Phase Approxima- 
tion (QRPA) and the nuclear shell model (NSM). They are in some sense 
complementary; QRPA uses a larger set of orbits, but truncates heavily the 
included configurations, while NSM can include only a rather small set of 
orbits but includes essentially all possible configurations. NSM also can be 
tested in a considerable detail by comparing to the nuclear spectroscopy 
data; in QRPA such comparisons are much more limited. 

For the 2v decay one can relate the various factors entering the cal- 
culations to other observables (f3 strength functions, cross sections of the 
charge-exchange reactions, etc.), accessible to the experiment. The consis- 
tency of the evaluation can be tested in that way. Of course, as pointed out 
above (see Fig ll.9|) the nuclear matrix elements for this mode are known 
anyway. Both methods are capable of describing the 2is matrix elements, at 
least qualitatively. These quantities, when expressed in natural units based 
on the sum rules, are very small. Hence their description depends on small 
components of the nuclear wave functions and is therefore challenging. In 
QRPA the agreement is achieved if the effective proton-neutron interaction 
coupling constant (usually called g pp ) is slightly (by ~ 10 - 20 %) adjusted. 

The theoretical description for the more interesting Ov mode cannot 
use any known nuclear observables, since there are no observables directly 
related to the M 0u . It is therefore much less clear how to properly esti- 
mate the uncertainty associated with the calculated values of M 0l/ , and to 
judge their accuracy. Since the calculations using QRPA are much simpler, 
an overwhelming majority of the published calculations uses that method. 
There are suggestions to use the spread of these published values of M 0v 
as a measure of uncertainty. 37 Following this, one would conclude that the 
uncertainty is quite large, a factor of three or as much as five. But that 
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way of assigning the uncertainty is questionable. Using all or most of the 
published values of M 0v means that one includes calculations of uneven 
quality. Some of them were devoted to the tests of various approximations, 
and concluded that they are not applicable. Some insist that other data, 
like the M 2v , are correctly reproduced, other do not pay any attention to 
such test. Also, different forms of the transition operator O 0v are used, in 
particular some works include approximately the effect of the short range 
nucleon-nucleon repulsion, while others neglect it. 

In contrast, in Ref. 38 an assesment of uncertainties in the matrix el- 
ements M 0v inherent in the QRPA was made, and it was concluded that 
with a consistent treatment the uncertainties are much less, perhaps only 
about 30% (see Fig |1.13[) . That calculation uses the known 2v matrix ele- 
ments in order to adjust the interaction constant mentioned above. There 
is a lively debate in the nuclear structure theory community, beyond the 
scope of these lectures, about this conclusion. 
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Fig. 1.13. Nuclear matrix elements and their variance for the indicated approximations 
(see Ref. 38 ). 



It is of interest also to compare the resulting matrix elements of Rodin 
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Table 1.1. Comparison of the calculated 
nuclear matrix elements M 0v using the 
QRPA method 38 and the NSM. 39 In the 
last column the NSM values are reduced, 
divided by 1.3, to account approximately 
for the effects of the induced nucleon cur- 
rents. 



Nucleus 


QRPA 


NSM 


NSM/1.3 


w Ge 


2.3-2.4 


2.35 


1.80 


82 Se 


1.9-2.1 


2.26 


1.74 


96 Zr 


0.3-0.4 






100 Mo 


1.1-1.2 






116 Cd 


1.2-1.4 






130 Te 


1.3 


2.13 


1.64 


136 Xe 


0.6-1.0 


1.77 


1.36 



et al. 38 based on QRPA and its generalizations, and those of the avail- 
able most recent NSM evaluation. 39 Note that the operators used in NSM 
evaluation do not include the induced nucleon currents that in QRPA re- 
duce the matrix element by about 30%. The QRPA 38 and NSM 39 M°" are 
compared in Table 11.11 In the last column the NSM matrix elements are 
reduced by 30% to approximately account for the missing terms in the op- 
erator, and to make the comparison more meaningful. With this reduction, 
it seems that QRPA results are a bit larger in the lighter nuclei and a bit 
smaller in the heavier ones than the NSM results, but basically within the 
30% uncertainty estimate. Once the NSM calculations for the intermediate 
mass nuclei 96 Zr, 100 Mo and 116 Cd become available, one can make a more 
meaningful comparison of the two methods. 

When comparing the results shown in Table 1X7X1 as well as the results of 
other calculations (e.g. Refs. 40,41 ) with Fig. 11.91 it is important to notice 
a qualitative difference in the behaviour of the 2v and Ois matrix elements 
when going from one nucleus to another one. For 2v the matrix elements 
change rapidly, but for the Ov the variation is much more gentle ( 96 Zr is 
a notable exception, at least for QRPA). That feature, common to most 
calculations, if verified, would help tremendously in comparing the results 
or constraints from one nucleus to another one. 

Once the nuclear matrix elements are fixed (by choosing your favorite set 
of results), they can be combined with the phase space factors (a complete 
list is available, e.g. in the monograph 42 ) to obtain a half-life prediction for 
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any value of the effective mass (mpp). It turns out that for a fixed (m^) 
the half-lives of different candidate nuclei do not differ very much from each 
other (not more than by factors ~ 3 or so) and, for example, the boundary 
between the degenerate and inverted hierarchy mass regions corresponds to 
half-lives ~ 10 27 years. Thus, the next generation of experiments, discussed 
above, should reach this region using several candidate nuclei, making the 
corresponding conclusions less nuclear model dependent. 

1.5. Neutrino magnetic moment and the distinction between 
Dirac and Majorana neutrinos 

Neutrino mass and magnetic moments are intimately related. In the or- 
thodox Standard Model neutrinos have a vanishing mass and magnetic mo- 
ments vanish as well. However, in the minimally extended SM containing 
gauge-singlet right-handed neutrinos the magnetic moment fi v is nonvan- 
ishing and proportional to the neutrino mass, but unobservably small, 43 

3eG -F o m-19 m " l-x m\ 

^ = V^ mV=3Xl ° ^TeV- (L19) 

Here hb is the electron Bohr magneton, traditionally used as unit also 
for the neutrino magnetic moments. An experimental observation of a 
magnetic moment larger than that given in eq. (|1.19[) would be an uneqivocal 
indication of physics beyond the minimally extended Standard Model. 

Laboratory searches for neutrino magnetic moments are typically based 
on the obsevation of the v — e scattering. Nonvanishing \x v will be recog- 
nizable only if the corresponding electromagnetic scattering cross section 
is at least comparable to the well understood weak interaction cross sec- 
tion. The magnitude of \x v (diagonal in flavor or transitional) which can be 
probed in this way is then given by 



1/2 

V^T^10- 10 f— ) , (1.20) 



Kl = G F m e 1 n _ 10 ( T 

Ms \/27i 



where T is the electron recoil kinetic energy. Considering realistic values 
of T, it would be difficult to reach sensitivities below ~ 10~ n /j.B- Present 
limits are about an order of magnitude larger than that. 

Limits on /i„ can also be obtained from bounds on the unobserved energy 
loss in astrophysical objects. For sufficiently large /i„ the rate of plasmon 
decay into the vv pairs would conflict with such bounds. Since plasmons 
can also decay weakly into the vv pairs , the sensitivity of this probe is 
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again limited by the size of the weak rate, leading to 

\tA = ^LHuir , (1.21) 
Mb V27ra 

where ujp is the plasmon frequency. Since usually (hiop) 2 <C m e T that 
limit is stronger than that given in eq. (|1.20|) . Current limits on /i„ based 
on such considerations are ~ 10~ 12 fiB- 

The interest in /i„ and its relation to neutrino mass dates from ~1990 
when it was suggested that the chlorine data 4 on solar neutrinos show an 
anticorrelation between the neutrino flux and the solar activity character- 
ized by the number of sunspots. A possible explanation was suggested in 
Ref. 44 where it was proposed that a magnetic moment fi u ~ 10~( 10-11 Vb 
would cause a precession in solar magnetic field of the neutrinos emitted ini- 
tially as left-handed v e into unobservable right-handed ones. Even though 
later analyses showed that the correlation with solar acivity does not exist, 
the possibility of a relatively large [i v accompanied by a small mass m v was 
widely discussed and various models accomplishing that were suggested. 

If a magnetic moment is generated by physics beyond the Standard 
Model (SM) at an energy scale A, we can generically express its value as 

^~X' (1 ' 22) 

where e is the electric charge and G contains a combination of coupling 
constants and loop factors. Removing the photon from the diagram gives 
a contribution to the neutrino mass of order 

m u ~ GA. (1.23) 

We thus arrive at the relationship 

m„~^^ ~ tttt^ [A(TeV)] 2 eV, (1.24) 

which implies that it is difficult to simultaneously reconcile a small neutrino 
mass and a large magnetic moment. 

This naive restriction given in eq. (|1.24p can be overcome via a careful 
choice for the new physics, e.g., by requiring certain additional symme- 
tries. 45-50 Note, however, that these symmetries are typically broken by 
Standard Model interactions. For Dirac neutrinos such symmetry (under 
which the left-handed neutrino and antineutrino v and v c transform as a 
doublet) is violated by SM gauge interactions. For Majorana neutrinos 
analogous symmetries are not broken by SM gauge interactions, but are 
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instead violated by SM Yukawa interactions, provided that the charged 
lepton masses are generated via the standard mechanism through Yukawa 
couplings to the SM Higgs boson. This suggests that the relation between 
\x v and m v is different for Dirac and Majorana neutrinos. This distinction 
can be, at least in principle, exploited experimentally, as shown below. 

Earlier, I have quoted the Ref. 17 (see Fig |1.4|) to stress that observation 
of the 0v/3f3 decay would necessarily imply the existence of a novanishing 
neutrino Majorana mass. Analogous considerations can be applied in this 
case. By calculating neutrino magnetic moment contributions to ra„ gen- 
erated by SM radiative corrections, one may obtain in this way general, 
"naturalness" upper limits on the size of neutrino magnetic moments by 
exploiting the experimental upper limits on the neutrino mass. 

In the case of Dirac neutrinos, a magnetic moment term will generically 
induce a radiative correction to the neutrino mass of order 51 



a A 2 \i v fj, v 



167rm e /!B 3x10 5 fiB 



A(TeV)] eV. (1.25) 



Taking A ~ 1 TeV and m„ < 0.3 eV, we obtain the limit /i„ < 10~ 15 /is 
(and a more stringent one for larger A), which is several orders of magnitude 
more constraining than current experimental upper limits on \x v . 

The case of Majorana neutrinos is more subtle, due to the relative flavor 
symmetries of m„ and /i„ respectively. For Majorana neutrinos the transi- 
tion magnetic moments [fJ-u] a (the only possible ones) are antisymmetric 
in the flavor indices {a,/3}, while the mass terms [m„] a ^ are symmetric. 
These different flavor symmetries play an important role in the limits, and 
are the origin of the difference between the magnetic moment constraints 
for Dirac and Majorana neutrinos. 

It has been shown in Ref. 52 that the constraints on Majorana neutrinos 
are significantly weaker than those for Dirac neutrinos, 51 as the different 
flavor symmetries of m v and /i„ lead to a mass term which is suppressed 
only by charged lepton masses. This conclusion was reached by considering 
one-loop mixing of the magnetic moment and mass operators generated 
by Standard Model interactions. The authors of Ref. 52 found that if a 
magnetic moment arises through a coupling of the neutrinos to the neutral 
component of the SU(2)l gauge boson, the constraints for fi Te and \i Ttl 
are comparable to present experiment limits, while the constraint on [i efi is 
significantly weaker. Thus, the analysis of Ref. 52 lead to a bound for the 
transition magnetic moment of Majorana neutrinos that is less stringent 
than present experimental limits. 
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Even more generally it was shown in Ref. 53 that two-loop matching 
of mass and magnetic moment operators implies stronger constraints than 
those obtained in 52 if the scale of the new physics A > 10 TeV. Moreover, 
these constraints apply to a magnetic moment generated by either the hy- 
percharge or SU(2)l gauge boson. In arriving at these conclusions, the 
most general set of operators that contribute at lowest order to the mass 
and magnetic moments of Majorana neutrinos was constructed, and model 
independent constraints which link the two were obtained. Thus the results 
of Ref. 53 imply completely model independent naturalness bound that - 
for A > 100 TeV - is stronger than present experimental limits (even for 
the weakest constrained element /i e/J ). On the other hand, for sufficiently 
low values of the scale A the known small values of the neutrino masses 
do not constrain the magnitude of the transition magnetic moment \i v for 
Majorana neutrinos more than the present experimental limits. Thus, if 
these conditions are fulfilled, the discovery of [i v might be forthcoming any 
day. 

The above result means that an experimental discovery of a magnetic 
moment near the present limits would signify that (i) neutrinos are Majo- 
rana fermions and (ii) new lepton number violating physics responsible for 
the generation of \i v arises at a scale A which is well below the see-saw scale. 
This would have, among other things, implications for the mechanism of 
the neutrinoless double beta decay and lepton flavor violation as discussed 
above and in Ref. 22 

1.6. Summary 

In these lectures I discussed the status of double beta decay, its relation to 
the charge conjugation symmetry of neutrinos and to the problem of the 
lepton number conservation in general. I have shown that if one makes 
the minimum assumption that the light neutrinos familiar from the oscilla- 
tion experiments which are interacting by the left-handed weak current are 
Majorana particles, then the rate of the Qvfi(3 decay can be related to the 
absolute scale of the neutrino mass in a straightforward way. 

On the other hand, it is also possible that the 0v/3(3 decay is mediated 
by the exchange of heavy particles. I explained that if the corresponding 
mass scale of such hypothetical particles is ~ 1 TeV, the corresponding 
Ov decay rate could be comparable to the decay rate associated with the 
exchange of a light neutrino. I further argued that the study of the lepton 
flavor violation involving /i — > e conversion and [i — > e + 7 decay may be 
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used as a "diagnostic tool" that could help to decide which of the possible 
mechanisms of the Ov decay is dominant. 

Further, I have shown that the the range of the effective masses (mpp) 
can be roughly divided into three regions of interest, each corresponding to a 
different neutrino mass pattern. The region of (rnpp) > 0.1 eV corresponds 
to the degenerate mass pattern. Its exploration is well advanced, and one 
can rather confidently expect that it will be explored by several (3(3 decay 
experiments in the next 3-5 years. This region of neutrino masses (or most 
of it) is also accessible to studies using the ordinary (3 decay and/or the 
observational cosmology. Thus, if the nature is kind enough to choose this 
mass pattern, we will have a multiple ways of exploring it. 

The region of 0.01 < (mpp) < 0.1 eV is often called the "inverted 
mass hierarchy" region. In fact, both the inverted and the quasi-degenerate 
but normal mass ordcrings arc possible in this case, and experimentally 
indistinguishable. Realistic plans to explore this region using the 0u(3(3 
decay exist, but correspond to a longer time scale of about 10 years. They 
require much larger, <~ ton size (3(3 sources and correspondingly even more 
stringent background suppression. 

Finally, the region {mpp) < 0.01 eV corresponds to the normal hierarchy 
only. There arc no realistic proposals at present to explore this mass region 
experimentally. 

Intimately related to the extraction of (mpp) from the decay rates is 
the problem of nuclear matrix elements. At present, there is no consensus 
among the nuclear theorists about their correct values, and the correspond- 
ing uncertainty. I argued that the uncertainty is less than some suggest, 
and that the closeness of the Quasiparticle Random Phase Approximation 
(QRPA) and Shell Model (NSM) results are encouraging. But this is still 
a problem that requires further improvements. 

In the last part I discussed the neutrino magnetic moments. I have 
shown that using the Standard Model radiative correction one can calcu- 
late the contribution of the magnetic moment to the neutrino mass. That 
contribution, naturally, should not exceed the experimental upper limit on 
the neutrino mass. Using this procedure one can show that the magnetic 
moment of Dirac neutrinos cannot exceed about 10~ 15 /xs, which is several 
orders of magnitudes less than the current experimental limits on \i v . On 
the other hand, due to the different symmetries of the magnetic moment 
and mass matrices for Majorana neutrinos, the corresponding constraints 
are much less restrictive, and do not exceed the current limits. Thus, a 
discovery of \i v near the present experimental limit would indicate that 
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neutrinos are Majorana particles, and the corresponding new physics scale 
is well below the GUT scale. 
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